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The stationary n-body problem in general relativity 

Robert Beig^ 
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Abstract. In this talk I describe recent joint works with R.Schoen and with 
' G. Gibbons and R.Schoen which prove the non-existence of certain asymptot- 

^ l' ically flat, stationary solutions of the Einstein equations with more than one 

, body. The basic restriction is for example satisfied when spacetime has an 

■ isometry reversing the sign of the timeUke Killing vector and fixing a hyper- 

' surface in the space of KiUing trajectories which is disjoint from the bodies. I 

' also give a detailed treatment of the Newtonian situation. 
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1. Introduction 



'A configuration of several gravitating bodies can not be in stationary equilib- 
rium': this plausible expectation is in the context of GR still awaiting a rigorous 
^ \ proof. We will in this talk start by explaining the Newtonian situation. This can 

' be summarized by the statement that there are no static solutions with more than 

one body, and where there exists a plane separating these bodies. This statement is 
' well-known, but I have not seen a printed version. That some separation condition 

\ is needed in the statement of that result is clear from the fact that without it there 

' are in fact solutions for the static 2 body problem, with source an elastic solid, 

\ both in the Newtonian [BSc| and in the Einstein theory |AS| . for example with 

one small body put near a critical point of the potential of a large body, this critical 
point lying in a suitably hollow region near the large body. 
. ■ In the next section we consider a nonlinear generalization of the Newtonian theory 

rS \ due to Giulini. Here we are able to obtain a weaker result in the spirit of the general 

i relativistic result in jBGSj . Namely, we show there can not exist solutions where 

the potential has zero normal derivative on a plane lying in the vacuum region. 
This result, which is in the spirit of the GR result in jBGS] . in particular rules out 
several-body solutions which are reflection symmetric across a plane in the vacuum 
region. In the final section we describe one of the results in jBGSj . This states 
that there can be no stationary spacetime having a reflection symmetry which in- 
verts the orientation of the orbits of the timelikc Killing vector and maps a timelike 
hypersurface tangential to these orbits into itself. This corresponds to a situation 
where the spin-spin interaction between matter on the two sides of this timelikc 
hypersurface is attractive. 

In the case of two axially symmetric black coaxial black holes the problem treated 
here has also been studied [WEj , and an analogous result has been proved in |LT| . 
In the absence of the discrete symmetry we are imposing, i.e. where there is the a 
priori possibility for spin-spin forces to balance the gravitational attraction, a full 
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proof of nonexistence has been found very recently in |NH| . using methods from 
the theory of completely integrable systems. 

2. The Newtonian situation 

Let ^(a) (ck — 1, 2, ..N) be bounded open connected sets with smooth boundary 
i9r2(Q) in R"^ with Euclidean metric 6ij. The sets ft(a) play the role of support 
of bodies. Each body has associated with it a positive function P(q;), the mass 
density, and the Cauchy stress tensor, a symmetric tensor field 'y(a)ij- In a- theory 
of continuum mechanics such as ideal elasticity, the quantities P(^a)i'\a)ij a-^c given, 
not as functions on M^, but as functions of some further fields such as maps from 
R'^ into some 'material manifold' Ba and of the first derivatives of these maps. The 
basic equations then take the following form (A = S'^^didj) 

(2.1) AU = inGj2P{a)Xn,^,, (7 ^ at oo 

a 

with > and 

(2.2) a(^a),j = P(a)dtU in il^^^ , '^(a)j"j = 

The detailed structure of the matter equations (j2.2p - which form a quasilinear 
second-order system of equations, usually elliptic, for the underlying map with 
Neumann-type boundary condition - need not occupy us here. The important thing 
is that these equations imply 

(2.3) M^i^,U) # / CdrU dx = 0, 

where is any of the 6 Euclidean Killing vector, i.e. satisfies 9(i = 0. The 
meaning of (j2.3p is of course that the total force and the total torque acting on each 
body be zero. One easily deduces (|2.3p from (j2.2p by using integration by parts, the 
symmetry of the Cauchy stress together with the Killing equation and the boundary 
condition of vanishing normal stress in (|2.2|) . 

The next question is if perhaps Eq. ()2.1|) alone implies these conditions. The answer 
is 'yes' for = 1. Namely the 'Newtonian stress tensor' defined by 

(2.4) e,, = ^ [{dmdjU) - ls,,{dur] 

satisfies, by virtue of (|2.ip . the relation 

(2.5) 9^e,,-^p(„)Xo<„,5.?7. 

a 

Thus the gravitational equation (j2.ip together with the matter equations (|2.2p - 
including the boundary conditions - imply that 

(2.6) 9^(0.j-E'^(")»J^"(»)) = O 

a 

holds in the sense of distributions. 

Integrating (|2.5p against over R'^ we get zero on the l.h. side: integrate by parts 
and use the Killing equation. The boundary term gives zero, since dU = 0(|x|~^). 
Thus 

(2.7) EM„(e,C/) = 

a 

In fact, defining ?7(q,) to be the unique solution of 

(2.8) AUi^a) = 47rGp(a) xn^^-, , U(a) -> at oo 

we find, by applying the same argument first to U^oi) and then to + t^(^), that 

(2.9) M„(^,[/(^))-fM^(e, [/(„)) =0 



THE STATIONARY N-BODY PROBLEM IN GENERAL RELATIVITY 



3 



This says that the total force and torque exerted by the gravitational field of body 
1 on body 2 is minus that exerted by body 2 on body 1 - which of course is nothing 
but the law 'actio = reactio'. Ea. (j2.9p also entails, for a = /3, the statement that 
the self-force and self-torque due to the gravitational field of a body vanishes. 
We now observe that (|2.9|) still falls short of guaranteeing p.3|) when iV > 1, 
because, for iV = 2 say, we need, in order for p.3|) to hold, in addition to Mi(^, [/(i)), 
M2{i, [/(2)) and U(2)) + M2(^, U(i)) all being zero, that also Mi(^, U(2)) and 

^2(CjC^(i)) vanish separately. 

Now, finally, we show that t/(2)) is indeed non-zero when ^{i) and f2(2) are 

separated by a plane. More precisely, when f2(i) lies to the left and ri(2) to the right 
of some separating plane 5", we will show that 

(2.10) 9„C/(2)|e>0 

on every plane E parallel to and to the left of S and where n is the left normal. This, 
by (|2.3p . will clearly imply Mi(^, C/(2)) > 0, when ^ is taken to be a translational 
Killing vector pointing along the left-normal of E. Inequality (|2.10p is nothing but 
the statement that the gravitational force due to body 2 on any body separated 
from body 2 by a plane is attractive. To prove (j2.10p . we use an elegant argument 
due to R.Schoen. Consider, in region 2, i.e. to the right of E, the field U{x) = 
U{2}{x) — U{2){p ° x), where p is the reflection across E. The field U satisfies 

(2.11) A^/ 47rGp(„)XO(„) in region 2 

Furthermore U goes to zero at infinity and vanishes on E. Thus, by the maximum 
principle, it follows that U < 0. Then, by the Hopf lemma, dnU\s = 2(?„C/(2)|s > 0. 
We have thus proved the 



Theorem 2.1. There exists no solution to the equations \2.1\2.2\) with N > I and 
where there is a plane S separating the bodies. 

Let us note that all the results obtained here could have been proved using the 
explicit Green's function of the fiat Laplacian, namely that, for compactly supported 
p, the unique solution U of AU = AnGp with U decaying at infinity has the form 

(2.12) U{x) = -g[ dx' 



We have avoided this and tried to make the argument as conceptual as possible. It is 
clear for example that the above result remains to be true, when the Laplacian A is 
replaced by 'meson- type' operator A — with ^ = const. We have however made 
crucial use of the linearity of the Poisson equation, which allows to consider the 
concept of 'force on a body due to the gravitational field of some other body'. This 
is not any longer available in the following nonlinear modification of the Newtonian 
theory. 

3. The Giulini theory 

This is a nonlinear modification of the Newtonian theory where 1/c^ of the 
gravitational self-energy has an active gravitational mass associated with it (the 
constant c being the speed of light). The equations now are 

(3.1) Af/-^ ^e(„)xn(„) - 2^(af/)2, C/^Oatoo 

(3.2) e^^d^ (e^cr(a)jj) = ^ d^U in , ^la)i^^^n^^-^ = 
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with e^jj) > playing the role of 'rest energy plus internal energy density'. The 
equation can be written in the linear form A^* = e(a)Xf2(„) in terms 

u 

of 4* = 6 2^, but that will play no role for us. The potential U has the expansion 

(3.3) u = -^ + 0^{^) 

The constant m in (j3.3p is positive. To prove this, observe that 

(3.4) lim / eii^ d^U dS' 

AttG R^^ Js„ 

where Sr is a large sphere. Now use the Gauss theorem together with (|3.1I) on the 
r.h. side of 

Again, in this theory, there is a law 'actio — reactio' in the form 

(3.5) 5^[e*(e,,-5]a(„),,xn<„,)] = 

a 

Now let S* be a plane disjoint from all bodies. We can then integrate p.5|) against 
over one of the half-spaces bounded by S, where ^ is the translational Killing 
vector coinciding with the unit normal n of 5*. Using integration by parts together 
with and ((^^ there results 

(3.6) i / [(a„t/)2 - {d''U){dAU)] dS = 

where Oa is a partial derivative tangential to S. Suppose that dnU\s = 0. It then 
follows that S is an equipotential for U. But this is impossible due to p.3p and the 
fact that m ^ 0. Thus we have obtained the 



Theorem 3.1. If a solution of i3.1l37^} contains a plane S disjoint from the bodies 
on which dnU is zero, this has to be the trivial ('no-body'-) solution with U ~ 0. 

The hypothesis of the above theorem is fulfilled when there is a reflection across 
some plane S disjoint from the bodies which maps the set of bodies into itself. It 
is this theorem - which is weaker than that available in the 'pure' Newtonian case 
- for which there exists an analogue in GR, to which we now turn. 

4. The GR situation 

The result I am now explaining works for the stationary Einstein equations with 
sources compactly supported in space. The nature of these sources is completely 
irrelevant except that they should meet the requirements of the positive-energy 
theorem, i.e. obey the dominant-energy conditiorQ. Otherwise the field equations 
will only be used in the vacuum region. For simplicity we will restrict to 3 -I- 1 
dimensions, although an analogous results for n -t- 1 with n > 3 is also true. We 
consider stationary - not just static - case, to allow for spin-spin interactions. Thus 
we assume (M, ds^) has a timclike Killing vector with complete orbits. Then ds^ 
can be written as (in distinction to the previous section, we set c equal to 1) 

(4.1) ds^ = -e^^{dt + i/jidx'f -f hijdx'dx^ 

with h being a Riemannian metric on N, the quotient space under the action of 
£, = dt- We refer to U as the gravitational potential and define 1/2 the curvature 
of the Sagnac connection i}} by 

(4.2) Lj='^d%p, ujij = d[iipj] 



Actually, all we need is that the ADM energy in the asymptotic rest system of the Killing 
vector be nonzero. 
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the vacuum field equations turn out to be {Qij = TZjj — ^hijTZ) 
(4.3a) -87rG(e.y +fi,y) =0 

(4.3b) AhU + e*^ ujkiuj''' = 

(4.3c) {e-^^ uj,j) = 

In (j4.3b ) we are using the definitions 

(4.4) 87rGe,j ^ 2 [iD,U)iD,U) - ^h,,{DUf] 
and 

(4.5) SttG = 2 e^^ [ -w^fco;/' + i /i^- ujkiuj'''] 

(Using du! = 0, Eq. (|4.3b ) follows from Eq.'s (|4.3b .b).) We consider asymptotically 
flat solutions of (|4.3p . They can be shown |BSij to have the form 



(4.6a) h^j = Sij + 0°° 



1 



(4.6b) u = -^ + 0-'(\ 

r V 



(4.6c) ujrj = -G • "T + O 
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in suitable coordinates. Here a' = — . Furthermore the constants m and i;,- — 
L[ij-^ are respectively the ADM mass and the (mass-centered) spin tensor of the 
configuration. When m = 0, the ADM energy of the initial data set induced ont = 
const is also zero. Then, when the source satisfies the dominant energy condition, 
the positive energy theorem implies that spacetime is fla10. 

Theorem 4.1. Let {N,h) have a hypersurface S disjoint from the matter region, 
which is non-compact, closed and totally geodesic w.r. to the unrescaled metric 
e~'^^ hij . Suppose in addition that the pull-hack-to-S of lo is zero. Then spacetime is 
flat. When spacetime has more than one asymptotically flat end, all the statements 
above hold separately w.r. to any such end. 

We first have to comment on the surface S in the hypothesis of the theorem. 
The typical situation where our hypothesis on S will be satisfied is when there is a 
reflection isometry 4* of spacetime which maps ^ into — ^ (and thusprojects down 
to an isometry of (TV, htj) and {N, e~^^hij)) and leaves S invariantQ 
In |BSc| it is shown in detail that S has to have a finite number of ends in the 
asymptotically flat region and S approaches a plane there in a precise sense. Thus, 
for each end of S asymptotic coordinates {x^,x^) in (|4.6[) can be chosen so that 
each end of S is given by = 0. It then follows from the hypothesis on uj that 
LjiB in (|4.6b ) vanishes. There are then good physical reasons detailed in [BGSj to 
believe that this corresponds to a situation where the spin-spin interaction between 
the bodies on the two sides of S is attractive. 

For the proof of the theorem we contract (|4.3b ) with n^n^ , with n* being the unit 
normal to S in {N, h). The result is 



(4.7) 



I {DnUf - i {DaU){D^U) 



AU A 
e UJAn^ n 



"^One can also allow for the presence of horizons, see IGHPI . IBC| . 

■^In fact the presence of the above isometry also implies DnU to vanish on S - a property we 
do not require in the theorem. 
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Here Dj^ is the intrinsic derivative on S. In (|4.7p we have used that wab is zero. 
Now. from the Gauss equation, 

(4.8) 3? = -2g„„ + (trfc)2-tr(A:2), 

where 3? is the Ricci scalar of S and k its extrinsic curvature. Since S is totally 
geodesic w.r. to e~'^^h, we have that 

(4.9) kAB ^ QAB DnU 

where q is the intrinsic metric on 5* induced from h. Inserting (|4.9l4.7p into Eq. (|4.8p . 
the terms involving DnU cancel so that finally 

(4.10) Ji^2iDAU)iD'^U) + 2e''^ ujAnU^^n ■ 

In particular "R is non-negative. One now integrates (|4.10p over S. Using the version 
of the Gauss-Bonnet formula for surfaces with boundary and the asymptotically 
planar nature of S, one arrives at a contradiction except when 31 vanishes and 
DaU, whence J7, is zero on S. Thus by (j4.6p m is zero, whence spacetime is flat, 
and we are done. 

We remark that for the n + 1-dimensional version of this theorem one uses, instead 
of the Gauss-Bonnet formula, the positive energy theorem in dimension n — 1 |SC| . 
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